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Abstract 

We consider the evolution of the light elements (Li, Be and B) incorporating the effects of 
their production by both neutrino process and cosmic-ray nucleosynthesis. We test the viability 
of the neutrino process to resolve the long standing problem of the 11 B/ 10 B isotopic ratio which 
amounts to 4 at the time of the formation of the solar system. This hypothesis may be ultimately 
constrained by the B/Be ratio observed in halo stars. Though we are able to obtain a solar 
isotopic ratio n B/ 10 B ~ 4, the current paucity of data at low metallicity prevents us from making 
a definitive conclusion regarding the resolution of this problem. We show however, that neutrino 
process nucleosynthesis leads to a relatively model independent prediction that the B / Be elemental 
ratio is large (> 50) at low metallicities ([Fe/H] < —3.0), if Be is produced as a secondary element 
(as is the case in the conventional scenario of galactic cosmic-ray nucleosynthesis). 



1 Introduction 



One of the successes of the standard big bang model is the prediction of the primordial abundances 
of the light nuclei, D, 3 He, 4 He, and 7 Li. In the standard model, the abundances of these isotopes, 
which span ten orders of magnitude (by number), are in accordance with observations for a narrow 
range in the baryon density or equivalently a baryon to photon ratio of 2.8 — 4 x 1CT 10 (Walker et 
al. 1991) . 

The other light nuclei, 6 Li, 9 Be, 10 B, and 11 B are thought to be produced in the interstellar 
medium by cosmic ray spallation and fusion reactions of protons and alpha particles with C, N, 
O and He nuclei. Their produced abundances in the standard big bang nucleosynthesis model is 
four to five orders of magnitude below the observations (Kajino & Boyd, 1990; Thomas et al., 
1993a). Production by cosmic rays can by and large account for the observed solar abundances 
of these elements (Reeves et al. 1970; Meneguzzi et al. 1971). However, 7 Li is also partly pro- 
duced by galactic cosmic-ray nucleosynthesis (GCRN) providing a potential constraint on big bang 
nucleosynthesis (Olive & Schramm 1992). Because ratios of the isotopes produced in GCRN are 
largely model independent, observations of the Be and B abundances can be an important test 
in determining what fraction of 7 Li observed in warm halo subdwarfs is primordial. Conventional 
GCRN models for the production of these elements which assume time-independent spectra and 
follow the metallicity of the ISM were able to explain the limited amount of data available at that 
time (i.e. for [Fe/H] > —1.4) (Vangioni-Flam et al. 1990). However these simple models tend to 
underproduce Be and B at low metallicities. These elements have now recently been observed in 
low metallicity population II stars (Rebolo et al. 1988; Ryan et al. 1990; Gilmore, Edvardsson & 
Nissen 1991; Ryan et al. 1992; Gilmore et al. 1992; Duncan, Lambert & Lemke 1992; Rebolo et al. 
1993). More recent models which treat the early galaxy as a "closed" rather than a "leaky" box for 
GCR increase the predicted abundances of Be and B for low metallicities almost to observed levels 
(Prantzos, Casse, & Vangioni-Flam 1992, hereafter PCV). 

It should be noted that inhomogeneities present in the early universe may alter the predictions of 
standard big bang nucleosynthesis (Applegate, Hogan & Scherrer 1987; Alcock, Fuller & Mathews 
1987) and give rise to increased nucleosynthesis beyond A = 8 ( Boyd & Kajino 1989; Malaney 
and Fowler 1989). However it has been shown (Thomas et al 1993b) that even under extreme 
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conditions the abundances of 6 Li, Be, and B in inhomogeneous big bang nucleosynthesis remain 
two to three orders of magnitude below their observed values in halo dwarf stars. Furthermore when 
the constraints from the light elements are applied, there is little distinction between standard and 
inhomogeneous model results for elements with A < 12 (Kurki-Suonio et al. 1990; Thomas et al. 
1993b; Reeves 1993). 

While the abundances of Li, Be, and B predicted by the improved models of GCR production of 
PCV almost match observations, these models cannot solve the long-standing problem of the solar 
U B/ 10 B ratio. Most GCR models predict a ratio of around 2.5 (Meneguzzi et al. 1971; Meneguzzi 
and Reeves 1975; Walker et al. 1985; Abia and Canal 1988; Steigman and Walker 1992; Walker 
et al. 1993). The observed ratio is found to be closer to 4 (Shima 1963; see also, Cameron 1982, 
Anders and Grevesse 1989). Meneguzzi and Reeves (1975) suggest as a possible solution to this 
problem introducing a steeply decreasing spectral component at low energies (less than 80 MeV/n) 
in GCR. At these energies, the cross sections of 10 B and U B are somewhat different, leading to a 
larger isotopic ratio. PCV point out, however, that a + a — > 6,7 Li reactions were ignored in these 
studies. PCV find that steepening the GCR spectra of protons and alphas at low energies leads to 
a solar isotopic ratio of U B/ 10 B ~ 4, but results in Li production considerably exceeding observed 
values in halo dwarfs. 

PCV also explore a possible solution to this problem by assuming that the more fragile 10 B is 
completely destroyed when ejected by a star of any size, while all of the U B is allowed to survive. 
Most previous models assume that all of the Li, Be, and B is destroyed when ejected from stars. 
The maximum isotopic ratio obtained under these conditions is 2.9. They conclude that this can 
not account for the observed isotopic ratio. 

With the difficulty in producing the observed isotopic ratio of n B/ 10 B, it has been suggested 
that alternative astrophysical sites for the production of 11 B must be found. One such site is at 
the shock front of type II supernovae, as suggested by Dearborn et al. (1989): when the shock hits 
the hydrogen envelope, it burns the ambient 3 He and A He producing 7 Be. Some of the resulting 
7 Be combines with alpha particles to produce 11 C which decays to 11 B. The primary goal of that 
work was to explore an alternative site for the production of 7 Li to reach Pop I abundances. They 
noted however that significant n B production might also take place. To date, the U B produced by 
this process has not been included in chemical evolution models and an isotopic ratio of U B/ 10 B 
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including production of 11 B by this process has not been determined. 

A potentially more important source for 11 B production has been found to result from neutrino 
induced nucleosynthesis in type II supernovae (Woosley et al. 1990). The core collapse of a massive 
star into a neutron star creates a flux of neutrinos so great that in spite of the small cross sections 
involved, it may still induce substantial nucleosynthesis. It is found that considerable 11 B production 
can result as the flux of neutrinos passes through the He, C, and Si shells of the stellar envelope, 
primarily by neutrino spallation of 12 C. In addition, some synthesis of 1 Li and 10 B takes place by 
this process but the production rate seems quite low. 

In this paper, we include production of ll B, 7 Li, and 10 B predicted from neutrino nucleosynthesis 
in type II supernova events (Woosley et al. 1993; Timmes et al. 1993) in galactic chemical evolution 
models along with GCR spallation production. We investigate the isotopic ratio of U B/ 10 B resulting 
from these models as well as the observational consequences of the increased abundances of B and 
Li. As we will show, although it is possible to obtain a U B/ 10 B isotopic ratio ~ 4, the B/Be ratio 
is typically high, B/Be >, 20 and is only marginally consistent with present data for halo stars. 
Indeed, this primary source of 11 B predicts significantly higher B/Be ratio at low metallicity than 
in standard models of GCRN and can as such conclusively test this hypothesis. 

The importance of the //-process was previously considered by Malaney (1992) primarily with 
regard to the production of 9 Be, though a discussion of the effect on B and the B / Be ratio was 
included. Malaney claimed that 9 'Be is produced via the 7 Li(t,n) 9 Be reaction. This rate was 
included in the recent Timmes et al. (1993) calculations and no substantial amount of 9 Be was 
found to survive the supernova shock. 

2 Galactic Evolution from the z/-process 

Before integrating the //-process yields into a full-blown model of galactic chemical evolution, we 
make the following, more simplified evaluation. Because we include a new source of B without any 
accompanying Be, we must be sensitive to the potential danger of an excessive B/Be ratio. In fact, 
this ratio has the possibility of strongly constraining the //-process source for B. The Be abundance 
observed in halo dwarf stars shows an almost linear dependence with metallicity. The observed Be 
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abundance can be fit approximately by 



[Be] = 12 + \og(Be/H) ~ 1.7 ± 0.4 + [Fe/i7] 



(1) 



Let us consider first the very simple model used in Steigman and Walker (1992) and Walker et 
al. (1993). There, it was assumed that the Pop II abundances of C,N and O in ISM were simply 
related to the iron abundance as [C/H] = [N/H] = [Fe/H] and [O/Fe] = 0.5. Evolution of these 
elements was not included in the model. Furthermore, it was also assumed that all exposure times 
for GCRN processes were equal. The overall exposure time then can be fixed by fitting the Be 
observations. Because the B/Be ratio in this model is relatively fixed at B/Be ~ 14 (i.e. the ratio 
of the corresponding spallation cross-sections), we can predict the B abundance as a function of 
[Fe], log(B/H) ~ —9.2 + [Fe/H] which is in good agreement with the fit to the three observations 
of B when a linear relation is assumed, 



The basis for the additional source of 11 B are neutrino processes occurring during type II super- 
novae. Neutrinos, including also and u T , are copiously produced in the hot collapsed core during 
a supernova (see eg., Mayle, Wilson and Schramm 1987). Because of their higher temperature, 
and v T neutral current reactions are dominant. The inelastic scattering of these neutrinos leads 
to unstable excited states which decay by p,n or a emission. These processes were included in 
supernova nucleosynthesis calculations by Woosley et al. (1990) where it was found that significant 
amounts of 7 Li and 11 B and lesser amounts of 10 B were produced (as well as other heavier isotopes). 

An important aspect of the calculation of Woosley et al. (1990) is the full treatment of pre- 
and post-shock nucleosynthesis. Since the duration of the neutrino burst exceeds the time scale for 
the passage of the shock through the inner layers of the exploding star, //-process nucleosynthesis 
can continue after the passage of the shock. In the outer layers, however, the destruction of fragile 
isotopes is a significant effect and is, for example, responsible for the destruction of 9 Be. Given 
the nature of the calculation there is a fair amount of uncertainty regarding the yields obtained 
by this process. The main uncertainty comes from the poorly known neutrino spectrum; indeed, 
there are indications that the high energy part of it (the most efficient for nucleosynthesis) should 
be considerably supressed (Myra & Burrows 1990). In this work, we include the recent results of 



[B] = 12 + hg(B/H) ~ 2.7 ± 0.9 + [Fe/H] 



(2) 
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Timmes et al. (1993) regarding the production of 7 Li, 10 B, and 11 B over a wide range of initial 
stellar masses. However, for the purposes of normalization to a solar n B/ 10 B ratio we take the 
liberty of adjusting the neutrino-process yields by an overall factor. 

To estimate the effect of the //-process on the boron isotopic ratio, we first separately integrate 
the //-process yields over a simple initial mass function with slope parameter, x = 1.8, and a 
constant star formation rate. We find as expected an approximately linear relationship for 11 B 
versus metallicity, [B] = 2.55 ± [Fe/H]. When added to the GCRN component for boron, we 
find that the U B/ 10 B isotopic ratio is raised from the canonical GCRN value of ~ 2.5 to 4.5. 
Furthermore, the B/Be ratio is increased from 14 to ~ 22. When compared with the available 
data: [B] = 0.4 ± 0.2 (HD 19445); = -0.1 ± 0.2 (HD 140283); and = 1.7 ± 0.4 (HD 201891) from 
Duncan et al. (1992) and [Be] < -0.3 (HD 19945) from Ryan et al. (1990); [Be] = -1.25 ± 0.4 
(HD 140283) from Ryan et al. (1992); = -0.97±0.25 (HD 140283) from Gilmore et al. (1992); and 
= 0.4 ±0.4 (HD 201891) from Rebolo et al. (1988). Thus the B/Be ratio becomes (propagating the 
quoted errors): B/Be = 20 ±26 (HD 201891); B/Be = 14 ±14 (HD 140283) based on the beryllium 
measurement of Ryan et al. (1992) and B/Be = 7± 5 (HD 140283) based on the measurements of 
Gilmore et al. (1992); and B/Be > 5 (HD 19445)0, one sees that although the agreement with the 
data is worsened when the //-process yield is included, this explanation of the boron isotopic ratio 
can by no means be excluded. 

Notice that in this section we assume a GCR origin for Be, 10 B and (part of) 11 B, but we 
adopt the observed B vs. Fe and Be vs. Fe relationships which are almost linear. It is well 
known however, that such linearity characterizes the evolution of a primary element (like oxygen), 
whereas in conventional GCR models B and Be behave as secondaries (see Prantzos et al. 1993b). 
This discrepancy between theory and the observations is somewhat alleviated, but not completely 
removed, in PCV, where a more efficient production of B and Be in the early Galaxy is suggested 
(see next section). Obviously, if the //-process B is incorporated in such a model, the resulting 
B / Be ratio at low metallicity is expected to be even larger than previously estimated. In the next 
section we investigate this effect with a fully developed model of galactic chemical evolution. 



to date only an upper limit for Be is available for this star (Ryan et al. (1990) 
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3 LiBeB production by neutrinos and cosmic rays 



We start with a brief description of the PCV model, where an important correction has been 
incorporated. It is currently thought that GCR are accelerated mainly by supernovae shock waves 
and propagate in the Galaxy by diffusing on the irregularities of the galactic magnetic field and 
suffering losses by ionization, nuclear reactions and leakage, in the framework of the "leaky box" 
model (e.g. Cesarsky 1980). Notice that in this model leakage takes place perpendicularly to the 
galactic disk, which has today a confinement volume with thickness H ~0.5 kpc and a gas fraction 
&gas,o ~ 0.10 - 0.20 (Rana 1991). The corresponding escape length (a key parameter caracterizing 
the leaky box model and indicating the amount of matter that GCR "traverse" before being lost) 
is estimated to be A ej0 ~ 10 g/cm 2 , at least up to 3-4 GeV/nucleon. 

In order to explain the recent B and Be observations in halo stars, PCV suggest that during 
its early evolution the Galaxy was a "closed" rather than a "leaky" box to GCR. According to the 
PCV analysis, the escape length evolves as A e oc a gas H, i.e. it was much larger in the past due to 
the larger gas content and the larger dimensions of the young Milky Way. The corresponding GCR 
spectrum is found to be flatter than the current one, and the associated proton flux above ~100 
MeV is found to evolve roughly as F oc VX e . Using a simplified prescription for the early collapse 
phase of the Galaxy and a detailed chemical evolution model, the recent Be and B observations 
are well reproduced with this hypothesis on GCR propagation. Notice, however, that the obtained 
B vs. Fe and Be vs. Fe relationships are still non linear, i.e. B and Be are clearly produced as 
secondary elements in that model. 

Since that work, Malaney and Butler (1993) noticed that the loss of protons due to ir° production 
has been neglected in PCV. This term is insignificant today. Indeed, the nuclear destruction length 
for protons is A^ ~ 200 gr cm" 2 and does not affect the total loss length A" 1 = A^ 1 + A" 1 ~ Aj 1 . 
However, this factor might have been important in the earliest phases of the Galaxy, if GCRs were 
as efficiently confined as suggested in PCV. In this work we incorporate the nuclear destruction 
term in our model. As a result we obtain GCR fluxes lower by a factor of ~2 w.r.t. the previous 
work at large A values (see Fig. 1 in Prantzos and Casse 1993). In the framework of the galactic 
evolution model this concerns quite a short period (a few 10 7 years), in the very beginning of the 
Galaxy. The resulting Be and B abundances are slightly affected, but still consistent with the 
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observations, as can be seen in Fig. 1 and 2. Notice also that the nuclear destruction term is much 
more important for alpha particles, which have a corresponding length of ~ 25 gr cm ~ 2 only. 
For that reason, GCR alphas never develop spectra as flat as GCR protons, even in the case of a 
good confinement: they are simply "removed" by collisions with the ISM protons for values of A > 
25 gr cm -2 . In this modified scenario (i.e. maintaining the efficient GCR confinement in the early 
Galaxy, but incorporating the nuclear destruction term) the early Li abundance is maintained low 
and consistent with observations (see Fig.3) The reason, however, is not the flattening of the GCR 
alpha spectra (as in the original PCV), but the nuclear destruction of alphas at large A. 

Given the reasonable success of the modified PCV model at explaining the low metallicity 
observations of Be and B, we now include new isotopic yields of Woosley et al. (1993) and Timmes 
et al. (1993) of //-process produced 7 Li, 10 B, and U B. Our goal is to explain of the U B/ 10 B 
ratio which is observed to be ~ 4 in the solar system. The PCV model, like almost all models of 
GCRN, gives a ratio of about 2.5. The additional source of U B from the //-process increases this 
ratio. We have adjusted the overall //-process yield to give a solar system U B/ 10 B ratio of 4. To 
do so, we found that it was necessary to reduce the neutrino induced yields by a factor ~ 2 to 
avoid overproducing U B. A summary of our results is displayed in Table 1. The resulting u B/ l0 B 
ratio as a function of metallicity is shown in Figure 4. This result can be compared directly to the 
case without a //-process contribution (dotted line) and to the case without an enhanced cosmic-ray 
spectrum (dashed line). As one can see, enhancing the cosmic-ray spectrum has little effect on the 
11 B/ 10 B ratio. The large discontinuity is solely due to our treatment of the two (halo and disk) 
separate components of the galaxy. (See PCV for details of the model.) A one zone model and 
presumably a more comprehensive treatment of the two zones (including infall in the disk) would 
show a smoother relationship. 

With the ll B/ l0 B ratio fixed, we can now examine the compatibility of the model with the data 
on Li, Be, and B. In Figures 1-3, we show the evolution of the LiBeB isotopes as a function of 
[Fe/H]. In each of these figures, we compare the abundance of the LiBeB element with (solid) and 
without (dotted) //-process nucleosynthesis and without an enhanced cosmic-ray spectrum (dashed). 
The solar values for B and Be are the same as those adopted in PVC. They are somewhat smaller 
than the values given by Anders & Grevesse (1989) and are closer to the values adopted in Reeves 
and Meyer (1978) and Arnould & Forestini (1989). In Figure 1, we show the beryllium abundance 
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as a function of [Fe/H]. As discussed above, the beryllium is slighlty lower than the results of 
PCV due to the correction in the escape length noticed by Malaney and Butler (1993), but is never- 
theless still in acceptable agreement with the data. The //-process makes very little contribution 
to Be/H. As shown in Figure 2, the model appears acceptable with respect to the evolution of 
the boron abundance both with and without the //-process contribution. For both beryllium and 
boron, the benefit of the enhanced cosmic-ray spectrum is clear. Finally we see that the lithium 
abundance shown in Figure 3 is also quite compatible with the data. Here, we have taken a 
primordial abundance [Li] ~ 2.0 (Olive and Schramm 1992) and have included both the GCRN 
and //-process contributions, neither of which is excessive at low metallicities. (Notice that since no 
stellar source of Li is included in the model, the predictions for the disk (—1 < [Fe/H] < 0) are 
lower that the observations.) 

Despite the apparent success regarding the observed abundances of the LiBeB elements, the 
B/Be ratio is more difficult to reproduce. In Figure 5, we show the B/Be ratio as a function of 
[Fe/H]. Though at [Fe/H] ~ — 1 to —2, the ratio is only slightly higher than standard GCRN 
predictions and consistent with the data for HD 201891, at lower [Fe/H] < —2 there is a clear 
departure from more standard model results. The B/Be becomes very large (B/Be > 50), for 
[Fe/H] <, — 3. This is clearly due to the primary nature of 11 B produced mostly by the //-process 
in the halo vs. the secondary nature of Be, produced only by GCRN. As it stands, the model is only 
marginally consistent with the data for HD 140283 (at the 2-a level). With regard to HD 140283, we 
note the following: 1) A recent observation (Molaro et al. 1993) has failed to identify beryllium in 
this star using a potentially more reliable part of the line spectrum and indicate that a lower value 
for Be/H and hence higher value for B/Be is likely; 2) The oxygen abundance in this star seems 
perhaps anomalously high. Several measurements indicate a value [O/Fe] ~ 0.8 — 0.9 (Bessel and 
Norris 1987; Molaro et al. 1993). Such a large oxygen enhancement may also be accompanied by a 
beryllium enhancement due to cosmic-ray spallation (the boron enhancement would be masked by 
the //-process contribution). Thus this star seems particularly incapable of excluding the high value 
for B/Be we would expect at [Fe/H] ~ —2.6 due to the //-process enhancement of 11 B. Future 
measurements of B and Be at low metallicity have the clear potential at confirming or rejecting 
the //-process contribution as a solution to the U B/ W B isotopic ratio. 
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4 Conclusion 



In summary, we have tested the hypothesis that a primary source of 11 B due to //-process nucleosyn- 
thesis in type II supernovae can explain the observed solar 11 B/ 10 B isotopic ratio. We note at this 
point that the observed value 11 B/ 10 B ~ 4 is due to a relatively old measurement (Shima 1963). 
Given the attention paid to this number, a new measurement of this ratio is certainly warranted. 
Though we are able to produce the oberved isotopic ratio of 4 as well as a reasonable evolution for 
the LiBeB elements, we expect a significantly higher B/Be ratio than in more standard models 
of cosmic-ray nucleosynthesis and galactic chemical evolution. Due to the primary nature of the 
neutrino-process source for U B relative to the secondary source for Be, we have found that the ratio 
B/Be > 50 when [Fe/H] < —3. Indeed, we think that measurements of the B/Be ratio in halo 
stars can help clarify the importance of neutrino-process nucleosynthesis (and ultimately confirm 
or invalidate it). 
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Table 1. Comparison between observed values at the birth of the solar system 

at t — 12.5 Gyr and the model 



Observed Data 



Model Result 




1 
1 
1 
1 



1 (normalization) 



1.03 
0.92 



1.08 



U B/ 10 B 
[O/H] 
[Fe/H] 
Z/Z Q 



4.05 ±0.1 





1 



4.11 (//-process yield adjusted) 



0.13 
0.011 
0.94 
0.19 



a (present value) 



0.10 - 0.20 
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Figure Captions 



Figure 1: The Be abundance as a function of [Fe/H] with (solid) and without (dotted) v- 
process nucleosynthesis included. The two curves coincide since 9 Be is not produced 
by the z/-process. Also shown is the Be abundance in an unenhanced GCRN model 
(dashed). The theoretical curves are compared to available data. (See text for 
references on observational data.) 

Figure 2: The B abundance as in Figure 1. 

Figure 3: The Li abundance as in Figure 1. Data are from Spite & Spite (1991). Only the 

upper envelope of the data is shown. 
Figure 4: The U B/ 10 B ratio as in Figure 1. 
Figure 5: The B/Be ratio as in Figure 1. 
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